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Abstract 

We present a detailed study of the derivation of the Hubbard model parameters for /c-(ET)2Cu2(CN)3 in the framework of ab initio 
Density Functional Theory. We show that calculations with different (i) wavefunction basis, (ii) exchange correlation functional 
and (iii) tight-binding models provide a reliable benchmark for the parameter values. We compare our results with available 
extended Hiickel molecular orbital calculations and discuss its implications forthe description of the properties of /(•-(ET)2Cu2(CN)3. 
D The electronic properties of ^•-(ET)2Cu(SCN)2 are also briefly discussed. 
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The search for possible candidates of a quantum spin liquid 
state has been one of the central issues in condensed matter 
physics since it was proposed several decades ago |1]. In the 
last years /c-(ET)2Cu2(CN)3 has attracted much attention due to 
its unique realization of a pressure-induced Mott transition from 
a spin liquid state [EllslSl to a metallic or even superconducting 
state at low temperatures. 

A large amount of theoretical investigations on the one-band 
half-filled Hubbard model on the anisotropic triangular lattice 
have been carried out in order to understand the mechanism 
of the phase transition. Irrespective of the current debate on 
the phase diagram, such as presence or absence of spin-liquid 
or superconducting state, derived from this simplified model 
by various methods, like, path integral renormalization group 
approach (5], variational Monte Carlo method |0, 01, cluster 
extension of dynamical mean-field theory |f^, or exact diago- 
nalization [9] a question of fundamental importance is whether 
the existing model parameters obtained from extended Hiickel 
molecular orbital calculations ifioll are reliable or not, as they 
are the basis of all the mentioned theoretical work. 

Recently we showed by full potential ab initio Density Func- 
tional Theory (DFT) calculations that the model parameters 
for /c-(ET)2Cu2(CN)3 have to be qualitatively revised from 
^7f=1.06 ^ to f7f=0.83+0.08 [11], which casts doubts on 
the validity of the description of the system as coupled quasi- 
one-dimensional chains. Similar results were obtained by the 
authors of Ref. |12]. In the present paper, we discuss the im- 
portance of including in the results for t' ft the effects of con- 
sidering different wavefunction bases as well as exchange cor- 
relation functionals within DFT. Furthermore, we show that 
alternative fitting procedures provide slightly different t' jt ra- 
tios. All these considerations lead to a statement of the t' It 
ratios within a benchmark. Finally we also discuss the effects 
of longer ranged hoppings. These considerations are important 
for all a:-(ET)2X charge transfer salts and we present here also 




Figure 1: Structure of the k BEDT-TTF arrangement for the example of k- 
(BEDT-TTF)2X with X=Cu2(CN)3. t and f' form a frustrated triangular lattice. 



details of the bandstructure calculations and parameter determi- 
nation for X=Cu(SCN)2. 

2. Crystal structure and bandstructure calculations 

The crystal structure of /c-(ET)2Cu2(CN)3 consists of alter- 
nate stackings of pairs of FT - BEDT-TTF molecule donors 
(see Fig. [U separated by Cu2(CN)3 acceptor layers. Each 
BEDT-TTF molecule consists of 10 carbon, 8 sulfur and 8 hy- 
drogen atoms linked by covalent bonds. A common feature 
of /f-(ET)2Cu2(CN)3 and other related BEDT-TTF based salts 
is the presence of conformational disorder. The two terminal 
ethylene groups in the BEDT-TTF molecule are twisted out of 
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Figure 2: Band structure and density of states of a:-(BEDT-TTF)2X with 
X=Cu2(CN)3, calculated with FPLAPW. The contribution of Cu to the DOS 
is highlighted. The Cu 3d states are fully occupied, indicating that Cu is in a 
nonmagnetic Cu*' oxidation state. 



Figure 3: Comparison of low energy band structure between different exchange 
and correlation potentials for /f-(BEDT-TTF)2X with X=Cu2(CN)3 . The FPLO 
basis set was used. For GGA, two Cu 3d bands are entangled with the bonding 
bands of the BEDT-TTF dimer molecules; for LDA, they are clearly separate 
(starting below about -0.6 eV) and therefore not shown. 



the molecular plane, and therefore this molecule has either an 
eclipsed or staggered conformation. The salt exhibits disorder 
at room temperature due to the random occurence of both con- 
formations. The /c-(ET)2Cu2(CN)3 salt crystallizes in the mon- 
oclinic space group P2\lc. The inversion centre symmetry in 
f 2] /c is guaranteed by the equal occupation of the two confor- 
mations. At low temperatures the staggered conformation dom- 
inates over the eclipsed conformation, therefore the staggered 
conformation was used in the calculations. The symmetry is 
then reduced to Pc with 120 atoms in the unit cell. 

Moreover, since the hydrogen positions were not avail- 
able from experiment [13] we introduced them in the struc- 
ture and relaxed their positions by performing Car-Parrinello 
projector-augmented wave molecular dynamics (CPMD) calcu- 
lations [15|,[16t]. For these calculations we considered a (4 
X 4 X 4) k-mesh in the Brillouin zone and plane wave cutoffs 
of 60 and 240 Ry for the wave function and the charge density, 
respectively. 

The electronic structure was analyzed by considering two 
different all-electron codes, namely the the Full Potential Lo- 
cal Orbital method (FPLO version 8.50-32) IH El] and the 
Full Potential Linearized Augmented Plane Wave (FPLAPW) 
method as implemented in Wien2k ifioll . We also performed 
comparative calculations with different exchange correlation 
functionals, the Local Density Approximation (LDA) in the 
parameterization of Perdew and Wang 1I20I1 and the Gener- 
alized Gradient Approximation (GGA) in the Perdew-Burke- 
Ernzerhof parameterization [21]. The convergency of the band 
structure with respect to the Brillouin zone integrations was 
checked with a series of calculations with up to 350 k-points 
in the full Brillouin zone. In the FPLAPW calculations a plane 
wave cut-off Rmt x ^max - 3.37 was used and the muffin- 
tin-radii for Cu, N, C, H, and S were chosen as 2.37, 1.06, 
1.16/1.06, 0.92, and 2.00 a.u., respectively. 

Despite the complexity of the crystal structure, the electronic 
structure of /^-(BEDT-TTF)2Cu2(CN)3 in the vicinity of the 



Fermi energy is surprisingly simple. Fig. |2] shows the band 
structure and density of states for /c-(BEDT-TTF)2Cu2(CN)3. 
Only two narrow, half filled bands with a width of about 0.4 eV 
cross the Fermi energy. Those bands are derived almost entirely 
from S 3p and C 2p states of the BEDT-TTF molecule and well 
separated from the remaining lower occupied and higher un- 
occupied bands. About 0.3 to 0.6 eV below the Fermi energy 
two further BEDT-TTF derived bands with mainly S 3p and C 
2p character are visible. This is in an overall good agreement 
with the results of extended Hiickel calculations [10], which 
find two anti-bonding bands of the highest occupied molecular 
orbital (HOMO) of the BEDT-TTF molecule intersecting the 
Fermi energy and place the corresponding bonding bands about 
0.4 eV below the Fermi energy. 

The Cu 3d states form a narrow band complex about 0.4 
to 2.1 eV below the Fermi energy. The Cu 3d shell is com- 
pletely filled, in agreement with the experimental reports [1(J 
that the Cu atoms are in the Cu^^ state with 10 d electrons. A 
comparison of the band structure of /c-(ET)2Cu2(CN)3 near the 
Fermi surface between the GGA and LDA functional is shown 
in Fig. [3] The two highest Cu 3d bands intersect the two bond- 
ing BEDT-TTF derived bands in GGA, whereas the LDA cal- 
culations place the Cu 3d bands further below the Fermi energy, 
so that the four BEDT-TTF derived bands near the Fermi energy 
are separated from the remaining bands within LDA. The differ- 
ence in the position of the Cu 3d level between LDA and GGA 
amounts to about 0.2 eV, resulting in a rigid band shift with 
only a minute effect on the shape of the BEDT-TTF bands near 
the Fermi energy and the model parameters derived from them 
(see below). A similar difference between LDA and GGA de- 
rived band structures is also found in related /(•-(BEDT-TTF)2-X 
compounds ifTlll . There are also minor differences in the band 
structures calculated with the two different electronic structure 
codes (Fig.Hl). While the two BEDT-TTF derived bands cross- 
ing the Fermi energy are almost exactly on top of each other, 
the lower occupied bands are shifted by about 50 me V between 
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Figure 4: Comparison of low energy band structure between FPLAPW and 
FPLO basis sets for k-(BEDT-TTF)2X with X=Cu2(CN)3 . The GGA exchange 
and correlation potential was used. 



0.2 



> 

0) 



0.2 



-0.4 



-0.6 



M 



Figure 5: Low energy bandstructure of k-(BEDT-TTF)2X with X=Cu(SCN)2 
(symbols) and tight binding fit (lines). Note he lifting of the degeneracy at the 
Z point. 



the FPLO and FPLAPW calculations. A possible reason for 
this difference could be the size of the basis set used in the cal- 
culations. However, we could not analyze the origin of these 
smaller differences in more detail, since the complexity of the 
crystal structure does not allow for extensive numerical checks. 

In addition to the Cu 3d bands there are further BEDT-TTF 
derived bands in the energy region between -2.5 and -0.75 eV 
giving rise to the complex band structure shown in Fig. |2] 
Seperated by a pseudo gap of about 0.5 eV a complex manifold 
of BEDT-TTF and anion layer derived bands is visible start- 
ing from 3 eV below the Fermi energy. The higher unoccupied 
bands around 2 eV above the Fermi level have predominatly S 
3p and C 2p character 

As comparison to the bandstructure of /<--(ET)2Cu2(CN)3, 
we present in Fig. |5] the low energy bandstructure for k- 
(ET)2Cu(SCN)2. While the overall features of the two bond- 
ing and antibondig bands are very similar to those of k- 
(ET)2Cu2(CN)3, /f-(ET)2Cu(SCN)2 shows a splitting of the 
bonding bands at the Z point which is not present in the band- 
structure of /(■-(ET)2Cu2(CN)3. 




Figure 6: Lattice of BEDT-TTF molecules (grey spheres) for /(■-(BEDT-TTF)2X 
with X=Cu(SCN)2 . Dimer bonds (ti ) are shown in black, f4 = t' in red. t has to 
be calculated from an average over ?2 (light green), (purple), fs (dark green) 
and ?6 (pink). The view is along the a direction, and the unit cell is marked. 



The lifting of the degeneracy at Z is due to the slightly dis- 
torted arrangement of FT molecules in /(•-(ET)2Cu(SCN)2 as 
shown in Fig. |6] This feature will also have implications on 
the determination of the tight-binding parameters as described 
below. 

3. Model Hamiltonian 

The half-filled Hubbard model 



H^y t(c^.c. +}i.c.) + y t'(c^. c. -hH.c.) 
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on the triangular lattice has been proposed as the suit- 
able Hamiltonian to describe the behavior of k-(BEDT- 
TTF)2Cu2(CN)3. f and f', as shown in Fig.[T]denote the hopping 
integrals between nearest neighbors < /, j > and next nearest 
neighbors [/, j] BEDT-TTF molecule dimers, respectively. Ex- 
tensions of this Hamiltonian include also longer ranged hopping 
integrals, like f", as well as nearest neighbor Coulomb integrals. 
In the following we shall concentrate on the calculation of the 
hopping matrix elements. 

In order to obtain the hopping integrals f and t' out of 
the bandstructure calculations presented above, we have per- 
formed various tight-binding studies. One option, which we 
denote dimer model, consists of defining a BEDT-TTF dimer 
of molecules as a lattice site. The corresponding tight-binding 
model is then fitted to the two antibonding bands in Figs. |3] 
andlH and the values of f, f', t" and longer-ranged hopping pa- 
rameters are directly obtained. The second option, or molecule 
model, considers the molecules in Fig.[T]as lattice sites and the 
corresponding tight-binding model is fitted to the four BEDT- 
TTF derived bands shown in Figs. |3] and |4] The parameters 
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f(meV) f' (meV) f It method 
Considering molecules as sites 
~^ 39 080 FPLO(GGA) 

50 42 0.83 FPLAPW(GGA) 

52 41 0.78 FPLO(LDA) 

Considering dimers as sites 
~50 45 090 FPLO(GGA) 

52 47 0.91 FPLAPW(GGA) 

53 47 0.88 FPLO(LDA) 

Table 1 : Hopping parameters t, t' and their ratio tjt' of /f-(ET)2Cu2 (CN)3 at 
ambient pressure, calculated with different basis sets (FPLO and FPLAPW) 
and functionals (LDA and GGA). 



extracted from this calculation are hopping integrals between 
molecules. In order to relate the molecule-molecule hopping 
parameters with the dimer-dimer hopping parameters (f, f , f" ) 
we consider the geometrical considerations given in Ref. m. 
This procedure was done for all bandstructure calculations dis- 
cussed in the previous section. In Table[T]we present the results 
for t and t' with the various approaches. 

While the differences between the sets of results are small, 
they are significant enough to be relevant for the determination 
of the f'/f ratio and define a margin of error of the DFT results. 
What is important about the results is that within the margin of 
error, the DFT calculations provide a clear prediction for the 
f'/f < 1 for /c-(ET)2Cu2(CN)3 in disagreement with extended 
Hiickel calculations that estimate t' jt > 1. Recent analysis of 
angular resolved magnetoresistance oscillations |22] predict a 
value of t' It for /(■-(ET)2Cu2(CN)3 of ~ 0.93 in very good agree- 
ment with our results. 

In the previous discussion we concentrated on the nearest and 
next nearest neighbor hopping integrals f, t' . It should be noted 
though that there are nonzero longer ranged hopping terms. In 
particular t" in the dimer model has a value of 7 meV it" jt - 
0.14), which is nonnegligible. 

In the tight-binding fit for /(■-(ET)2Cu(SCN)2 the lifting of the 
degeneracy at the Z point can only be captured if the f2, f3 and 
fs, fg of Fig.|6]are differentiated even though ti ~ and t^ ^ t^. 

Finally we would like to note that while the half-filled Hub- 
bard model on an anisotropic triangular lattice with t and f hop- 
ping parameters is the minimal model that we can consider for 
the description of the system, a few extensions may be impor- 
tant like: (i) inclusion of longer ranged hopping terms, (ii) in- 
clusion of longer ranged Coulomb interaction terms and (iii) 
consideration of the system at the level of molecules in terms 
of the quarter-filled Hubbard model ll23ll . 



culations. This analysis shows that the model parameters 
for /(•-(ET)2Cu2(CN)3 have to be qualitatively revised from 
t'lt=\M ^ to f7f=0.83+0.08. Implications for the proper- 
ties of /(■-(ET)2Cu2(CN)3 and possible extensions of the models 
are discussed. 
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4. Conclusions 

In summary we have presented a detailed study of the 
derivation of Hubbard model parameters for /('-(ET)2Cu2(CN)3 
in the framework of Density Functional Theory. Compari- 
son of the results obtained with different exchange and cor- 
relation functionals, high precision electronic structure codes 
as well as different tight binding models allows us to estab- 
lish well defined bounds for the t' ft ratio used in model cal- 
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